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Abstract 

Variable annuities (VA) are popular insurance products. VAs provides 
the insured with a guaranteed accumulation rate on their premium at 
maturity. In addition, the insured may receive extra benefit if returns 
of underlying funds are high enough. Here we consider a special case of 
VA with high-water mark feature and Guaranteed Minimum payment re- 
set. In Black-Scholes model for underlying fund we derive explicit pricing 
formula for this type of contract. The value of VA contracts depends 
on the time between observation dates. Corrections due to this effect 
are calculated and compared with Monte-Carlo results. Good agreement 
between analytical formula and numerical calculations of VA values is 
demonstrated. 

1 Introduction 

Variable annuity (VA) contract provides the policyholder with a guaranteed 
minimum return and offers extra benefit if returns of underlying funds are high 
enough. VAs with high water mark feature are similar to lookback option: 
they choose the maximal fund value over the entire term of an annuity for 
calculating payoff at maturity. Within Black-Scholes framework [1], these type 
of options were considered in [2], [3], [4]. The effects of stochastic interest rates 
and mortality were explored in [5] . 

In this paper, we focus on valuation of VA contract value where guaran- 
teed minimum payoff is calculated from set of fund values over the entire term 
of contract. High- water mark values are calculated at predetermined observa- 
tion dates. Usually, frequencies of these reset dates are fixed and they can be 
monthly, quarterly or yearly. Because of this feature, some of maximal fund 
values can be missed and it reduces value of the contract. This effect needs to 
be taken into account to price correctly VA contract. 
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2 Contract Payoff Definition 



We consider the following lookback contract, which pay-off at expiration time 
T can be presented in the following form 

Payoff = V(T) = NMax{Hw{T), (1 + 1 At)V{t - At)) (1) 

where AT is a notional 1 ; V(0) = 1; 



Hw(t) = Max(S(t)/S(0),Hw(t - At)); 

Hw(0) = 1; 



(2) 



is a High Watermark; S(t) is price of underlying at time t; At is time interval 
between observation dates; 7 is a rate of guaranteed minimum payoff. 

In continuous limit At equation (1) can be presented in the following 
form 



V(T) 



Max 
e lT Max 



Sit). 
5(0) 



e^ T -V\t G [0,T] 



m 

5(0) 



e-^lt G [0,T] 



(3) 



From eq.(3) we can see that the probability distribution of the VA contract 
value is distribution of extreme (maximal) values of underlying 



S(t) = 5(i)e~ 7t 



(4) 



So, we can use distribution of extreme (maximal) values to calculate present 
value of this VA contract. 



3 Extreme Value Distribution in Black-Scholes 
Model 

Consider the following dynamics for underlying S(t) 

d\n(S(t)) = fxdt + vdW(t) (5) 

where 

1 2 

fi = r-y--v (6) 

is a risk-neutral drift; r and y are short term interest rate and dividend yield; v 
is volatility; W{t) is a standard Brownian motion 

< dW(ti)dW(t 2 ) >= S(ti - t 2 )dt (7) 

1 Below we assume that Af = 1 and 5(0) = 1 
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Notice, that to model distribution of underlying S(i) (4), all we need to do, is 
just to redefine dividend yield as 

y^y + i (8) 

Probability distribution function of returns 

x(t) = HS(t)/S(0)) (9) 

in Black-Scholcs world (5) satisfies the following Fokker-Plank equation 

dP(x,t) _ v 2 d 2 P{x,t) dP(x,t) 

dt ~ 2 dx 2 M dx { U) 

with initial boundary condition 

P(x,0) = 5(x) (11) 

This solution is 

^o(x,t) = -=i=e- 1 ^T i (12) 

To find probability distribution function of maximal values x(t) for time 
interval < t < T it is convenient to calculate probability of event that x(t) < h 
for all t in the time interval [0,T]. This problem is equivalent to problem with 
absorbing wall at x = h and therefore we need to find solution which satisfies 
eq.(10) for all < t < T and x < h and to impose the following boundary 
conditions: 

P(x = h,t) = 0. (13) 

Solution of eq.(13) with boundary conditions (13) can be found from the fol- 
lowing anzatz: 

P(x,t,h) = P (x,t) + AP (x - x ,t) (14) 

where A and xq are constants and xq > h. 

Eqs. (13,14) lead us to the following equation 

P(h,t,h)= e ^r\\-Ae~~^e )= 0. (15) 

and we obtain that 

A = e 2 $\ x = 2h. (16) 

Then, probability of the event, that return x does not cross the border x = h 
by time t = T can be calculated 



P(x <h)= ( P(x,h,T)da 

J —OO 
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where 



2 f z 2 
erf(z) = —= \ e~ x dx 
V 71 " Jo 

erfc(z) = 1 - erf(z) (18) 

are error functions. 

It gives us the following formula for probability distribution of maximal 
values 

P(h,fi,v,T) = ^P(x<h) = 

2 _ (h-Tuf (J, 2t H (h + Tji\ 

= e 2v^t -e ~ eric , 19 

4 High Watermark Liability Value 

In the limit At — > present value of VA contract (3) can be calculated from 
maximal value distribution (19) as 

V{S,X,T) = e h ~ r)T Max [S(t) e -^\t e [0,T},X] = 

= e^- r}T X + e^- r ^ T / (Se h - X)P(h, fi, v, T)dh = 

Jln(X/S) 

°"> T (Xe- rt + 2 Call(S, X, v, r, y + 7 , t)+ 



e 



v 2 + 2/i 



M(l + ar/( ^±^Mj7a )),^ + 

l^\ 1+ ^ vV/c f^±ln(A7^)^ 
V \/2^t / 



(20) 



where /j, = r — y — ^ ; X = Sh > S(0) is a strike; 5jy is maximal observed 
price of underlying from issue date to current time (t = 0). Call{S, X, v, r, y, t) 
is a Black-Scholes formula for European Call Option price. Here, to calculate 
present value of that contract we take into account discount factor e~ rT . 

High Watermark contract value depends on size of time interval between 
observation dates (At > 0). Value of contract with At > is lower than in the 
limit At — > just because selected dates of resets (observation dates) can miss 
the actual extremal values of underlying stochastic process. 

To take into account finite size of time steps, let us consider Brownian motion 
(5) on time interval T with N equal time intervals between observation dates: 

0, t\, t2, ■ ■ ■ , t n , . . . , T; 

t n = n dt; dt = T/N. (21) 
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Then, the average maximal value near the point U is 

(Max(e x{t) \t e [U - At/2, t t + At/2], 

\ * / h 

= e x[u) (l + v\ 2- + 0(dt)) (22) 



It means that in the limit of small time steps At — > 0, distributions of extreme 
values can be presented in the following form 

P(h, (J., v, T, At) ~ P(h, n, v.T, At = 0) (23) 

where 

h~h + v\T^-. (24) 
Then, the integral in (20) can be transformed to the following form 

/>oo 

/ (Se h - X)P(h,n,v,T,At)dh~ 

poo 

~ / (Se h - X)P(h + e,fi,v,T, At = 0)dh = 

J\n(X/S) 
poo 

= e- e (Se h - X)P(h,fi,v,T,At = Q)dh (25) 

J\n(X/S) 

where e = v^foM] X = Xe e . 

y i" 

And finally, from (20) we obtain the following formula for VA contract value 

p OO 

XeP' r)T + e~ rT / (Se h - X)P{h, ft, v, T, At)dh ~ 



XeP^ T + e 



pOO 

/ (S"e' 1 - X)P(h,n,v,T,At = Q)dh (26) 

J\n(X/S) 



where X = 5'#e 7th ; t^ is a time from the observed underlying value to current 
time which gives the highest guaranteed payoff at time t = 0. 
At issue date, the initial contract value is 

poo 

S e (7-r)T + e -rTg / ( e h _ yp( h) M) W) Tj A ^ d/l _ 

JO) 
/>oo 

~ S*e (7 ~ r)T + e~ e S / (e /l - l)P(/i,^,w,T,At = 0)dh (27) 
Jo 

Fig.(l) demostrates a good agreement between analytical formula (27) and nu- 
merical calculations of VA values. 



5 



1.7 -i 



Eq.(20) 

Eq.(27) 



W "I 1 1 v 

10 100 1000 10000 

Number of Observation Dates 



Figure 1: VA Contract Value vs Number of Observation Dates. T = 10, r = 5%, 

7 = 8%, « = 10% 

5 Conclusion 

In this paper we consider VA contract which guaranteed minimum payoff is 
calculated from a set of maximal fund values over the entire term of contract. 
In the limit of small time intervals between observation dates (At — > 0) we 
derived explicit pricing formula for this type of VA contracts. Finite time in- 
terval corrections are calculated and compared with Monte-Carlo results. Good 
agreement between analytical formula and numerical calculations of VA values 
is demonstrated. 
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